home *** CD-ROM | disk | FTP | other *** search
/ Libris Britannia 4 / science library(b).zip / science library(b) / DJGPP / LGP250S1.ZIP / src / libgplus.5 / libgplus / tests / trationa.cc < prev    next >
C/C++ Source or Header  |  1993-06-06  |  5KB  |  192 lines
  1. /*
  2.  a test file for Rational class
  3.  */
  4.  
  5. #include <Rational.h>
  6.  
  7. #include <assert.h>
  8.  
  9. #define tassert(ex) {if ((ex)) cerr << #ex << "\n"; \
  10.                        else _assert(#ex, __FILE__,__LINE__); }
  11.  
  12.  
  13. void identtest(Rational& a, Rational& b, Rational& c)
  14. {
  15.   Rational one(1, 1);
  16.  
  17.   assert(-(-a) == a);
  18.   assert((a + b) ==  (b + a));
  19.   assert((a + (-b)) ==  (a - b));
  20.   assert((a * b) ==  (b * a));
  21.   assert((a * (-b)) == -(a * b));
  22.   assert((a / (-b)) == -(a / b));
  23.   assert((a / b) == (a * (one / b)));
  24.   assert((a / b) == (one / (b / a)));
  25.   assert((a - b) ==  -(b - a));
  26.   assert((a + (b + c)) == ((a + b) + c));
  27.   assert((a * (b * c)) == ((a * b) * c));
  28.   assert((a * (b + c)) == ((a * b) + (a * c)));
  29.   assert(((a - b) + b) == a);
  30.   assert(((a + b) - b) == a);
  31.   assert(((a * b) / b) == a);
  32.   assert(((a / b) * b) == a);
  33.  
  34.   Rational x = a;
  35.   x *= b;
  36.   assert(x == (a * b));
  37.   x += c;
  38.   assert(x == ((a * b) + c));
  39.   x -= a;
  40.   assert(x == (((a * b) + c) - a));
  41.   x /= b;
  42.   assert(x == ((((a * b) + c) - a) / b));
  43.  
  44.   assert(x.OK());
  45. }
  46.  
  47.  
  48.  
  49. void simpletest()
  50. {
  51.   Rational one = 1;
  52.   assert(one.OK());
  53.   Rational third(1, 3);
  54.   assert(third.OK());
  55.   Rational half(1, 2);
  56.   assert(half.OK());
  57.  
  58.   Rational two(2);
  59.   Rational zero(0);
  60.   Rational r;
  61.   r = two+zero;
  62.  
  63.   cout << "one = " << one << "\n";
  64.   cout << "two = " << r << "\n";
  65.   cout << "third = " << third << "\n";
  66.   cout << "half = " << half << "\n";
  67.  
  68.   cout << "third + half = " << third + half << "\n";
  69.   cout << "third - half = " << third - half << "\n";
  70.   cout << "third * half = " << third * half << "\n";
  71.   cout << "third / half = " << third / half << "\n";
  72.  
  73.   Rational onePointTwo = 1.2;
  74.   cout << "onePointTwo = " << onePointTwo << "\n";
  75.   cout << "double(onePointTwo) = " << double(onePointTwo) << "\n";
  76.  
  77.   Rational a = one;
  78.   cout << "a = " << a << "\n";
  79.   assert(a.OK());
  80.   a += half;
  81.   cout << "a += half = " << a << "\n";
  82.   assert(a == Rational(3, 2));
  83.   a -= half;
  84.   cout << "a -= half = " << a << "\n";
  85.   assert(a == Rational(1));
  86.   a *= half;
  87.   cout << "a *= half = " << a << "\n";
  88.   assert(a == half);
  89.   a /= half;
  90.   cout << "a /= half = " << a << "\n";
  91.   assert(a == Rational(1));
  92.   assert(a.OK());
  93.  
  94.   identtest(one, one, one);
  95.   identtest(one, third, half);
  96.   identtest(third, half, one);
  97.   identtest(onePointTwo, half, a);
  98. }
  99.  
  100. void pitest()
  101. {
  102.   Rational half(1, 2);
  103.   Rational approxpi(355, 113);
  104.   assert(approxpi.OK());
  105.   cout << "approxpi = " << approxpi << "\n";
  106.   cout << "double(approxpi) = " << double(approxpi) << "\n";
  107.  
  108.   Rational rpi = Rational(M_PI);
  109.   cout << "rpi = Rational(PI) = " << rpi << "\n";
  110.   assert(rpi.OK());
  111.   cout << "double(rpi) = " << double(rpi) << "\n";
  112.  
  113.   cout << "approxpi + rpi = " << approxpi + rpi << "\n";
  114.   cout << "approxpi - rpi = " << approxpi - rpi << "\n";
  115.   cout << "approxpi * rpi = " << approxpi * rpi << "\n";
  116.   cout << "approxpi / rpi = " << approxpi / rpi << "\n";
  117.  
  118.   Rational negapproxpi = -approxpi;
  119.  
  120.   cout << "-approxpi = " << negapproxpi << "\n";
  121.   assert(sign(negapproxpi) < 0);
  122.   cout << "abs(negapproxpi) = " << abs(negapproxpi) << "\n";
  123.   assert(abs(negapproxpi) == approxpi);
  124.  
  125.   assert(approxpi != rpi);
  126.   assert(approxpi >= rpi);
  127.   assert(approxpi > rpi);
  128.   assert(!(approxpi == rpi));
  129.   assert(!(approxpi <= rpi));
  130.   assert(!(approxpi < rpi));
  131. #ifdef __GNUC__
  132.   assert((approxpi >? rpi) == approxpi);
  133.   assert((approxpi <? rpi) == rpi);
  134. #endif
  135.  
  136.   assert(floor(approxpi) == 3);
  137.   assert(ceil(approxpi) == 4);
  138.   assert(trunc(approxpi) == 3);
  139.   assert(round(approxpi) == 3);
  140.  
  141.   assert(floor(negapproxpi + half) == -3);
  142.   assert(ceil(negapproxpi + half) == -2);
  143.   assert(trunc(negapproxpi + half) == -2);
  144.   assert(round(negapproxpi + half) == -3);
  145.  
  146.   identtest(approxpi, rpi, negapproxpi);
  147.   identtest(rpi, approxpi, rpi);
  148.   identtest(negapproxpi, half, rpi);
  149. }
  150.  
  151.  
  152.  
  153. void IOtest()
  154. {
  155.   Rational a;
  156.   cout << "\nenter a Rational in form a/b or a: ";
  157.   cin >> a;
  158.   cout << "number = " << a << "\n";
  159.   assert(a.OK());
  160. }
  161.  
  162. // as a fct just to test Rational fcts
  163. Rational estimate_e(long n)
  164. {
  165.   Rational x = Rational(n + 1, n);
  166.   Rational e = pow(x, n);
  167.   return e;
  168. }
  169.  
  170. void etest(long n)
  171. {
  172.   cout << "approximating e as pow(1+1/n),n) for n =" << n << "\n";
  173.   Rational approxe = estimate_e(n);
  174.   assert(approxe.OK());
  175.   cout << "double(approxe) = " << double(approxe) << "\n";
  176.   cout << "log(approxe) = " << log(approxe) << "\n";
  177.   assert(log(approxe) <= 1.0);
  178.   cout << "approxe = " << approxe << "\n";
  179. }
  180.  
  181. int main()
  182. {
  183.   simpletest();
  184.   pitest();
  185.   IOtest();
  186.   etest(10);
  187.   etest(100);
  188.   etest(1000);
  189.   cout << "\nEnd of test\n";
  190.   return 0;
  191. }
  192.